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The adsorption dynamics of a colloidal particle at a fluid interface is studied theoretically and 
numerically, documenting distinctly different relaxation regimes. The adsorption of a perfectly 
smooth particle is characterized by a fast exponential relaxation to thermodynamic equilibrium 
where the interfacial free energy has a minimum. The short relaxation time is given by the ratio of 
viscous damping to capillary forces. Physical and/or chemical heterogeneities in a colloidal system, 
however, can result in multiple minima of the free energy giving rise to metastability. In the presence 
of metastable states we observe a crossover to a slow logarithmic relaxation reminiscent of physical 
aging in glassy systems. The long relaxation time is determined by the thermally-activated escape 
rate from metastable states. Analytical expressions derived in this work yield quantitative agreement 
with molecular dynamics simulations and recent experimental observations. This work provides new 
insights on the adsorption dynamics of colloidal particles at fluid interfaces. 
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The adsorption and binding of colloidal particles to 
fluid interfaces is relevant to numerous natural and in- 
dustrial processes. Novel technological applications in 
areas that range from materials science to renewable en- 
ergy [U [2] , and from food science to biomedicine [3| l4| 
demand advancements in the fundamental understanding 
of colloidal adsorption dynamics. Standard models based 
on continuum thermodynamics [l] |5] predict a monotonic 
relaxation to an equilibrium position where the contact 
angle with the interface is given by Young's law [6 . Ac- 
cording to standard models, this equilibrium position cor- 
responds to a stable state determined by the (global) 
minimum of the Helmholtz free energy. For micrometer- 
or nanometer-size particles, the energy decrease at the 
equilibrium state can be orders of magnitude larger than 
the thermal energy, and thus strong interfacial forces are 
expected to cause a spontaneous adsorption with rapid 
relaxation to equilibrium. 

Nevertheless, the dynamics of colloidal adsorption re- 
mains poorly understood for systems of great practical 
interest (e.g. functionalized colloidal particles). Funda- 
mental issues arise when the equilibrium contact angle is 
difficult to determine, e.g. due to the presence of physical 
and/or chemical heterogeneities [6 , 7 . Furthermore, even 
when the equilibrium contact angle can be properly de- 
termined it is frequently observed that the adsorption of a 
colloidal particle is neither fast nor spontaneous, and re- 
quires some form of external actuation (mechanical, ther- 
mal, or chemical) to be initiated |8HTQ]. Notably, recent 
experimental work [TO^ reported a slow logarithmic relax- 
ation to equilibrium after initiating the adsorption of a 
micrometer-size particle at a water-oil interface. The un- 
expected observation was attributed to nanoscale surface 
heterogeneities [10]. A logarithmic relaxation is reminis- 
cent of physical aging in glassy systems having complex 
free energy landscapes and metastability pT] . 



These phenomena suggest that the presence of en- 
ergy barriers associated with microscale heterogeneities 
should be considered in order to describe the adsorption 
dynamics of colloidal particles. In this Letter, we study 
the dynamics of adsorption in the presence of multiple 
metastable states caused by spatial fluctuations in the 
interfacial free energy. We propose a model that quanti- 
tatively describes recent experimental observations [lOj. 
Our derived analytical expressions and numerical sim- 
ulations reveal a non-trivial adsorption dynamics with 
crossovers from logarithmic to exponential relaxation. 

We begin our analysis with the equation of motion for 
a colloidal particle that straddles the interface between 
two fluids as illustrated in Fig. [TJa). Assuming that the 
particle undergoes Brownian motion, and neglecting the 
action of external fields, its vertical position z is governed 
by a Langevin equation 

mz = y/2k^r]{t) - + F{z); (1) 

here m is the particle mass, ksT is the thermal energy of 
the surrounding fluids, r]{t) is a zero- mean unit- variance 
Gaussian noise, ^ is the viscous friction coefficient, and 
F{z) = —dU/dz is the interfacial or capillary force de- 
termined by the interfacial free energy U . In the frame- 
work of continuum thermodynamics, the energy to form 
an interface is the product of the interfacial area and the 
corresponding surface tensions: 7 for the fluid-fluid inter- 
face; 7pi and 7^2 for the interfaces between the particle 
and each fluid phase. Hence, for a spherical particle of 
radius R with its center at vertical position z and a sharp 
fluid interface located at z = 0, the interfacial free energy 
can be cast as [5j 

Us{z) = \K{z - ZE? - C for \z\ < R, (2) 

where z^; is the equilibrium position, while K = 27r7 and 
C = 7Tj{R — zeY positive constants. According to 
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Eq.|2]the interfacial force is linear, F = —K{z — ze)- Em- 
ploying this linear force and neglecting small inertial ef- 
fects, the solution of Eq. [l] gives an average displacement 
< z >= ze + Az ex.p{—t/T]j)^ where Az — z(0) — ze is 
the distance from equilibrium at t = and Td = ^/i^ is 
the viscous decay time. For reference, we must note that 
the decay time to equilibrium is ~ 0.1/is for a one- 
micron radius particle adsorbed at a water-oil interface. 
Thus, the standard model predicts a fast exponential re- 
laxation to equilibrium for a perfectly smooth particle. 

A few comments are in order. For the idealized case of 
a perfectly smooth and spherical particle that straddles 
a flat interface, there is a one-to-one correspondence be- 
tween the particle position z and the observed contact an- 
gle d = acos(— 2;/i?) which is constant along the circular 
contact line located at z = 0. According to Young's law 
for the equilibrium contact angle, cos Oe = (7^2 — 7pi)/7, 
the particle will straddle the two fluids at an equilib- 
rium position Ze = —RcosOe when \jp2 — 7^1 1 < 7. 
Both Young's law and Eq. |2]for Us{z) apply to homoge- 
neous and regular contact lines. In order to model hetero- 
geneities that cause fluctuations of the contact line shape 
and its mean position, we will neglect the line tension 
term used in alternative models [121 US] and introduce 
spatial fluctuations in the interfacial free energy. 

We are interested in a simple analytical description 
of microscale physical/chemical heterogeneities, as illus- 
trated in Fig. [ijb), that can give rise to significant en- 
ergy barriers AU > ksT. For this purpose, we introduce 
spatial fluctuations of the interfacial free energy with the 
general form ^AU sin (A^ +(/)); here ^AU is the amplitude 
of the fluctuation, / = 27r/A is its wavelength, and <p is 
an arbitrary phase. For the sake of analytical tract abil- 
ity we will consider a single-mode fluctuation of small 
wavelength, / <C equilibrium contact angles near neu- 
tral wetting, 70° ^ Oe ^ 110°, and conditions where the 
particle center is close to the fluid interface, \z/R\ <C 1. 
Curvature effects, of order 0{{z/R)'^)^ can thus be ne- 
glected and the free energy is expressed as 

U{z) = Us{z) + ^AUsm{Xz + (j)) (3) 

while the capillary force on the particle is 

F{z) = -K{z - Ze) - ^XAU cos{Xz + (/)). (4) 

A simplified system with dynamics governed by Eq. [l] 
and free energy given by Eq. [s] is depicted in Fig. [ijc). 
The simplified system consists of a solid surface element 
with sinusoidal roughness (wavelength /, height /i, and 
width w) that is in contact with the fluid interface, and 
thus produces interfacial energy fluctuations of magni- 
tude AU = jhw as it undergoes Brownian motion. Just 
as if the surface element formed part of a much larger 
spherical particle, a linear force drives the system toward 
equilibrium z ^ ze] this force corresponds to the free 
energy contribution /J^ in Eq. [2] caused by the macro- 
scopic curvature of the particle. The decomposition of 




FIG. 1. Problem description and modeled system, (a) Col- 
loidal particle of radius R straddling a sharp fluid interface, 
(b) Particle with modeled microscale heterogeneities of size 
/ <C R. Interfacial energy fluctuations have wavelength / and 
amplitude AU = 'jhw associated with surface roughness of 
height h and width w ^ h. The system free energy is given by 
Eq. [3] (c) Simplified system for MD simulations with free en- 
ergy given by Eq. [S] A linear "spring" force F = KAz drives 
the relaxation to equilibrium at z = ze- (d) Trajectories z{t) 
from MD simulations (dotted lines) for different realizations 
{LH=Kl/2kBT=7A, AU/kBT=4, and Az=-4.5Lh)- Analyt- 
ical expressions are plotted for comparison (see legend). 



macroscale and microscale features makes the dynamics 
tractable via molecular dynamics (MD) simulation [14]. 
Simulated trajectories z{t) from individual MD realiza- 
tions shown in Fig. [ijd) exhibit several metastable states 
with a lifetime that increases as the equilibrium position 
Ze is approached. 

According to Eqs.[3]-[4] there are multiple minima in the 
free energy for \z — ze\ < ttAU /Kl, and sufficiently close 
to equilibrium the system must exhibit metastability. 
The basins of attraction of each metastable state are cen- 
tered at the local minima Zo = ZE-\-l{Ti—^—(j)/27r)-\-0{e) 
and each basin is bounded by two neighboring maxima 
at z± = ± ^/ + 0(e); here n is any integer and 
e = Kl\z — zeI/ttAU is a small parameter. A particle 
undergoing Brownian motion (see Fig. [TJd)) will transi- 
tion, or "hop" back and forth, between metastable states 
at a local rate 15 






d^U{z±) 







exp 



(5) 



predicted for \z — Zo\ < \L The energy barriers, AU± = 
U{z±) — U{zo)^ in the forward/backward direction de- 
termine the Arrhenius exponential factor. The prefac- 
tor employed in Eq. [5] is valid for overdamped systems 
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[151 ES] where ^ > ^m\d'^U{z±) / dz'^\ (here m is the 
particle mass). For \z — ze\ <C ti/SU/KI the motion 
is dominated by thermahy-activated hopping, and the 
ensemble- averaged speed < i > of the particle is deter- 
mined by an ordinary differential equation 



<z>= ^/(r+ -r_). 



(6) 



We solve analytically the equation of motion in Eq. |6] to 
obtain the average trajectory 



{z) = ze + Lh log 



1 + AgexpH/Tg) 
[l-Anex^i-t/TH) 



(7) 



which depends on three independent parameters: the 
characteristic hop length Lh = 2kBT/Kl; the trajec- 
tory "amplitude" Ah = tdinh{^Az /Lh) determined by 
the initial separation from equilibrium Az = z{0) — ze 
at time t = 0; and the characteristic hop time 



Th = T^ 
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1J_ 



(8) 



In Eq. Is] we introduced the ratio <l> = \AU\^/K be- 
tween the free energy curvature of the modeled sinu- 
soidal fluctuation and that of a smooth spherical par- 
ticle. The characteristic hop length, Lh^ and time, Tf/, 
in our model are thus determined by two independent pa- 
rameters, (i) the dimensionless energy barrier AU/kBT 
and (ii) the dimensionless wavelength 1/ ^ksT /^^ while 
the dynamics of relaxation also depends on the dimen- 
sionless initial condition Az/Lh- Sufficiently far from 
equilibrium when \z — ze\/Lh ^ 1, we have {z) = 
Ze + Ln^ogl^t/TH + e-xjp{Az/ Lh)]^ which is equivalent 
in form to the logarithmic expression recently employed 
in Ref. [10] to fit experimental observations by treating 
Lh and Th as adjustable parameters. It is noteworthy 
that Eq. [T] predicts a crossover to exponential relaxation, 
{z) = Ze -\- Azexp(— t/Tjc/), for \z — ze\/Lh <C 1 when 
the particle is very close to equilibrium. 

To verify our analytical predictions, derived under the 
assumption of a quasi-static process and (overdamped) 
Brownian motion [TSl [16 , we perform numerical simu- 
lations via standard MD techniques The simulated 
system, illustrated in Figjljc), has three atomic species: 
{i=l) fluid 1; {i = 2) fluid 2; {i = 3) the solid parti- 
cle and bottom wall. As in previous work our MD 
simulations employ generalized Lennard- Jones potentials 



V{r) = 4e[(r/cr) 
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Cij{r/a) where e is the inter- 



action energy, a is roughly the atomic radius, r is the 
distance between any two atoms, and Cij = Cji is the 
interaction coefficient between species (z, j = 1-3). In 
this work we set ca = 1 for self interactions, while in the 
case of cross interactions we set Ci2 = 0.5 for the fluids, 
Ci3 = C23 = 0.35 for the fluids and Brownian particle, 
and Ci3 = C23 = 0.8 for the fluids and stationary bottom 
wall. At a simulated constant temperature T = 3e//cB, 
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FIG. 2. Mean relaxation trajectories {z(t) — ze)- Hori- 
zontal dotted lines indicate the distance from equilibrium 
Azh = \tiAU/1 above which hopping dominates. Markers 
indicate numerical simulations via MD and LD (see legend). 
Solid lines indicate analytical predictions from Eq. [7| Dashed 
lines indicate exponential relaxations at the viscous decay 
time Td = ^/K and the hop time Th from Eq. 8] (see la- 
bels). The energy barrier is AU — AksT^ and the fluctuation 
wavelength is / = A.byJksT /'^ for all three initial conditions: 
(a)Az/Lij=-1.8; (b) Az/Lif=-4.6; and (c) Az/Lij=-7.3. 



maintained by a Nose-Hoover thermostat, and a mass 
density p = 0.8/cr^, the fluids are macroscopically im- 
miscible and the surface tension measured across a plane 
interface [19 is 7 lAkBT/cr^ . Solid surfaces exhibit 
neutral wetting {Oe = 90°) given the symmetry of fluid- 
solid interactions. All atoms have a unit mass and are 
initialized on a fee lattice (cf. Fig. [ijc)) with spatial 
spacing Ax = y^l/p. The solid particle and wall, carved 
from the fee lattice, are neutrally buoyant. The parti- 
cle has length Lt = 40Ax and a sinusoidal roughness 
with wavelength / = 6Ax, height h = 2Ax, and width 
w = 2/1, which in turn produces surface energy fluctu- 
ations of amplitude AU = ^yhw. In different numerical 
realizations, fluid atoms are initialized with random ve- 
locities and the particle is allowed to move in the vertical 
z direction after thermal equilibrium is attained. The 
hop length Lh = 2kBT/Kl in MD simulations is mod- 
ified by adjusting the "spring" stiffness K. In addition, 
we perform Langevin dynamics (LD) simulations which 
are equivalent to the numerical solution of Eq. [T] 

Ensemble-averaged trajectories {z{t)) are reported in 
Fig. [2] for three different displacement amplitudes \Az\ = 
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FIG. 3. Crossover from exponential to logarithmic relax- 
ation. Solid lines indicate logarithmic relaxation given by 
Eq. [71 Dashed lines correspond to an exponential decay at 
rate 1/Td. Dotted horizontal lines show the distance to equi- 
librium /\zh = ^TvAU/Kl near which the crossover occurs. 



\z{0) — ze\ — 2-7Lh' The simulated case corresponds to 
an energy barrier AU = AksT and a fluctuation wave- 
length / = 4.5y^/c^T/7. Very close to equilibrium for 
\Az\/Lh < 2 (cf. Fig. [2|a)), the relaxation is exponen- 
tial at the slower rate 1/Th predicted by Eq.[7| Far from 
equilibrium (cf. Figs. [2|b-c)) where \z — ze\ > — 
\i\A\J jl^ the relaxation is exponential at the "fast" rate 
I/Td predicted by Eqs. [l]-[2] for a stable system; To is 
numerically computed using large values of K for which 
there is no metastability. Closer to equilibrium where 
\z — ze\ < ^zh^ numerical results are in close agree- 
ment with Eq.[7| valid for thermally- activated transitions 
between met ast able states. 

Additional simulations for larger energy barriers 
AU = 10-20/cbT, different fluctuation wavelengths 
/ = 4-16y^TbT77, and a larger separation from equi- 
librium |A2:| = 32 = 20-82Lh are reported in 
Fig. [3] The resulting relaxation times Th in the logarith- 
mic regime are extremely long and these simulations were 
only feasible via LD. At a distance Azh = ^ttAU/KI 
we observe the crossover from a fast exponential relax- 
ation, driven by surface energy minimization, to a slow 
logarithmic relaxation predicted by Eq. ^ driven by the 
thermally-activated escape from metastable states. The 
crossover to a slow logarithmic relaxation is delayed, or 
even prevented, when increasing the wavelength / of the 
fluctuation, cf. Figs[3];^a-b), or decreasing the energy bar- 
rier AU, cf. Fig.[3|;a) and Fig.[3];c). 

It is useful to examine how features of the en- 
ergy fluctuations and physical properties of the me- 
dia determine the relaxation time Th in the loga- 
rithmic regime. According to Eq. [S] the dimension- 



less hop time Th/Td is a function of the dimension- 
less amplitude and wavelength of the energy fluctuation 
Th/Td = f{AU/kBT, l/^ksTh). Given the functional 
form of Eq. [Sj it is convenient to analyze the function 
\og(TH/TD) = P-Q, where P = AU/ksT + /ksT 
is the dominant contribution. One flnds that not only 
increasing the energy barrier At/, but also increasing the 
wavelength / causes an exponential increase in the re- 
laxation time Extremely large relaxation times Th 
produced by large fluctuation wavelengths / > ^Jk^Tj^, 
indicate that the particle will be jammed, i.e. prevented 
from reaching equilibrium, as soon as the hopping motion 
begins. However, when the wavelength / of the hetero- 
geneity is large the particle can get much closer to equilib- 
rium before jamming because the hopping motion begins 
at a much later stage when \z — ze\ — ^zh ^ AU/Kl. 
Therefore, a logarithmic relaxation on experimentally ac- 
cessible time scales can only be observed when the energy 
fluctuation wavelength / < ^JkB^^J^ is smaller than the 
mean-square displacement of the fluid interface. 

In order to discuss our flndings let us consider a phys- 
ical system of practical importance, a one-micron radius 
and nearly spherical particle adsorbed at a water-oil in- 
terface for which 7 = 0.04N/m and ksT = 4 x 10"^^J 
at room temperature. The contact line perimeter is then 
p :^ 6 X 10~^m. We further consider a moderate en- 
ergy barrier AU = jAA = 30/cbT, caused by a sur- 
face heterogeneity of area A A 3 nm^. If the defect 
is well localized only a small portion of the contact line 
hops over the defect while most of its perimeter remains 
pinned. Hopping over localized defects, the contact line 
moves in steps of average length A/p :^ 0.5 x 10~^^m. 
This distance corresponds to an energy fluctuation wave- 
length / 0.002 y^/cbT/7 which results in a relaxation 
time Th — lO^s. The logarithmic regime will then begin 
at \z — ze\ ^ 5 X 10~^m, i.e. half a radius away from the 
expected equilibrium. This scenario seems to describe 
the experimental conditions studied in Ref. [10 . More- 
over, using a wavelength / 0.5-1 xlO~^^m and energy 
barriers AU 15-30/c^T, we flnd that Eq. [7|flts closely 
the experimental data reported in Ref. [10 . 

In conclusion, we have derived an analytical expres- 
sion for the adsorption dynamics of a colloidal particle in 
the case of thermally-activated relaxation to equilibrium. 
We found that moderate energy barriers, of order 10k bT, 
due to localized nanoscale heterogeneities can give rise 
to physical aging as the particle approaches the expected 
equilibrium position ze- Numerical simulations show a 
crossover from a fast exponential relaxation to a slow log- 
arithmic relaxation at a distance \z — ze\ ~ AU/jl from 
the expected equilibrium. These results demonstrate a 
nontrivial relaxation dynamics of adsorption of colloidal 
particles at a fluid interface. 
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